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Abstract 

We study the effects of the unstable gravitino on the big-bang nucleosynthesis. 
If the gravitino mass is smaller than 10 TeV, primordial gravitinos produced after 
the inflation are likely to decay after the big-bang nucleosynthesis starts, and the 
light element abundances may be significantly affected by the hadro- and photo- 
dissociation processes as well as by the p *^ n conversion process. We calculate the 
light element abundances and derived upper bound on the reheating temperature 
after the infiation. In our analysis, we calculate the decay parameters of the gravitino 
(i.e., lifetime and branching ratios) in detail. In addition, we performed a systematic 
study of the hadron spectrum produced by the gravitino decay, taking account of 
all the hadrons produced by the decay products of the gravitino (including the 
daughter superparticles). We discuss the model-dependence of the upper bound on 
the reheating temperature. 



1 Introduction 



Low-energy supersymmetry, which is one of the most prominent candidates of the physics 
beyond the standard model, may significantly affect the evolution of the universe. One 
reason is that, assuming i?-parity conservation, the lightest superparticle (LSP) is stable, 
which becomes a well-motivated candidate of the cold dark matter. Another famous rea- 
son, which is very important in the framework of local supersymmetry (i.e., supergravity), 
is that there exist various possible very weakly interacting particles in the particle con- 
tent. The most important example is the gravitino, which is the gauge field for the local 
supersymmetry. 

Since the gravitino is the superpartner of the graviton, its interaction is suppressed 
by inverse powers of the gravitational scale and hence its interaction is very weak. Even 
though the gravitino is very weakly interacting, however, it can be produced by the scat- 
tering processes of the standard-model particles (and their superpartners) in the thermal 
bath in the early universe. Once produced, the gravitino decays with very long lifetime. 
In particular, if the gravitino mass is smaller than ~ 20 TeV, the lifetime becomes longer 
than ~ 1 sec and hence the primordial gravitinos decay after the big-bang nucleosynthesis 
(BBN) starts. The (unstable) gravitino is expected to be relatively heavy, and its decay re- 
leases energetic particles which cause dissociation processes of the light elements generated 
by the standard BBN reactions. Since the standard BBN scenario more or less predicts 
light-element abundances consistent with the observations, those dissociation processes are 
harmful and, if the abundance of the primordial gravitinos are too much, resultant abun- 
dances of the light elements become inconsistent with the observations. Such argument 
provides significant constraints on the primordial abundance of the gravitino [1^.*^ 

In the inflationary scenario, which is also strongly suggested by the WMAP data |3], 
gravitinos are once diluted by the entropy production after the inflation. Even in this case, 
however, gravitinos are produced again by the scattering processes of the particles in the 
thermal bath. Since the total amount of the gravitinos produced by the scattering pro- 
cesses is approximately proportional to the reheating temperature Tr, the BBN scenario 
is too much affected to be consistent with the observations if the reheating temperature 
is high. In the past, the BBN constraints on the unstable gravitino have been intensively 
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Recently, Kawasaki and two of the present authors (K.K. and T.M.) have studied the 
constraints on the unstable long-lived particle from the BBN in detail P^IT^ : in particular, 
in this paper, effects of the hadrons produced by the decay of such unstable particles (as 
well as the effects of the photo-dissociation) were systematically studied,*^ and general 

Constraints on the case with the gravitino LSP have been considered in Refs. [21IS1- 
*^In order to relax the constraints, several scenarios have been studied. In Ref. USj, it was discussed that 
the modification of the expansion rate by the 5D efi^ect in the brane-world cosmology, which may reduce 
the abundance of the gravitino. In Ref. they studied possibilities of the dilution of the gravitino by 
the late-time entropy production due to the decaying moduli without newly producing many gravitinos. 
"^^^For old studies on the effects of the hadronic decay modes, see also |2()l I21j . 
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constraints on the primordial abundance of such unstable particles were presented. Then, 
the results were applied to the case of the unstable gravitino and the upper bound on the 
reheating temperature was obtained for several simple cases. 

In P^ll6j. however, several simplifications and assumptions are made for the properties 
of the gravitino. First of all, several very simple decay patterns of the gravitino were 
considered, which are applicable for very specific mass spectrum of the superparticles. 
In addition, for the hadronic branching ratio of the gravitino, only several typical values 
were used to obtain the constraint. Furthermore, for some cases (in particular, for the 
case where the gravitino dominantly decays into the gluon and the gluino), effects of 
the hadrons emitted from the superparticles (like the gluino) were neglected. Thus, it 
is desirable to perform more detailed and complete analysis of the upper bound on the 
reheating temperature with the unstable gravitino. 

In this paper, we study the effects of the unstable gravitinos on the BBN, paying 
particular attention to the properties of the gravitino. Compared to the previous works, 
decay processes of the gravitino (and the decay chains of the decay products including the 
superparticles) are precisely and systematically studied. As a result, energy spectra of the 
hadrons (in particular, proton, neutron, and pions) produced by the gravitino decay are 
studied in detail for various mass spectrum of the superparticles. 

The organization of this paper is as follows. In Section |2l we present the model we 
consider and summarize the important parameters for our analysis. In Section 01 detail 
of the decay processes of the gravitino is discussed. Some of the important issues in our 
analysis, which is the secondary decays of the daughter particles from the gravitino decay 
and their hadronization processes, are discussed in Section |3J Then, in Section outline 
of our calculation of the light-element abundances is discussed. Our main results are given 
in Section ini and Section [7| is devoted for conclusions and discussion. 

2 Model 

In this paper, we adopt the minimal particle content to derive the upper bound on the 
reheating temperature. Thus, the model we consider includes the particles in the minimal 
super symmetric standard model (MSSM) as well as the gravitino. These particles are 
listed in Table [TJ 

In order to precisely calculate the decay rate of the gravitino, it is necessary to ob- 
tain the mass eigenvalues and mixing parameters of the superparticles. Thus, we briefly 
summarize the relation between the gauge-eigenstate and mass-eigenstate bases here. 

We start with the neutralino sector. With the SU{2)l and U{1)y gaugino masses M2 
and Ml as well as the supersymmetric Higgs mass /i^, the mass matrix of the neutralinos 
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Particles Notation 



Gravitino ipfj, 
Neutralinos xl, xl, xl 

Charginos xf , xt 

Gluino g 

Squarks q = ml, mr, dh, ^r, 5l, 5r, cl, cr, 6i, 62, h, ^2 

Sleptons Cl, cr, /xl, /^r, n, f2, z>eL, u^^, z/^^ 

Photon 7 

Weak bosons Z, 

Gluon g 

Neutral CP even Higgses h, H 

CP odd Higgs A 

Charged Higgs 

Quarks q = u, d, s, c, b, t 

Leptons e, /x, r, Ve, J^fj., 

Table 1: List of particles in the mass-eigenstate bases. 



is given in the form*^ 



( Ml —mz sin 6w cos (3 mz sin 6y/ sin (3 \ 

M2 mz cos 6w cos [3 —mz cos 6^1^ sin [3 

-mz sin cos (3 mz cos 6*^/ cos /? — /^i? 

Y mz sin 6'vi/ sin /3 —mz cos ^vf sin /3 — /^h y 



(2.1) 



where mz is the Z-boson mass while 9w is the weak mixing angle, and tan (3 is the ratio 
of the vacuum expectation value of up- and down-type Higgs bosons. This mass matrix is 
diagonalized by a unitary matrix f/^o as 



U*oM^oUj = diag(m^o,m^o,m^o,m^o). 
In addition, for the chargino sector, the mass matrix is given by 



M2 V2mw cos /3 
\^mw sin /3 fj^n 



(2.2) 



(2.3) 



where mw is the Ty^-boson mass. We define the unitary matrices diagonalizing Ai^^ 

U*^+M^±U-} = diag(m^±,m^±). (2.4) 



and Uy^ 



^^We used the convention of 
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For the neutral Higgs boson, the gauge-eigenstates (i.e., the up-type Higgs and 
down-type Higgs and the mass-eigenstates are related by using the mixing angle a: 

(H\_ cosa sin a \ / Re(i/0) - \ , . 

\h) ~^\-smacosa) \Re{H^^) - V2 ) ' ^ ' 

where Vi and V2 are vacuum expectation values of and respectively. 

In addition, we have to consider the mixings in the squark and slepton mass matri- 
ces. For simplicity, we do not consider the generation mixing in the squark and slepton 
sector. We also neglect the left-right mixing in the first and second generation squarks 
and sleptons since such mixing is small in many class of models, in particular, in the 
minimal supergravity (mSUGRA) type models which we adopt in our analysis. For the 
third-generation squarks and sleptons, we take account of the effects of the left-right mix- 
ing. Such mixing is parameterized by unitary matrices C/j, C/^, and Uf which diagonalize 
the mass-squared matrices of the squarks and sleptons: 

U^M'fJj^ = dmg{m\,m\) : f = i,b,f. (2.6) 

In our analysis, wc consider the case where the gravitino is unstable and the LSP is 
contained in the MSSM particles. Since the charged or colored LSP is disfavored, we 
consider the case where the LSP is the lightest ncutralino Xi- Consequently, the gravitino 
is assumed to be heavier than Xv Of course, the gravitino may be heavier than other 
superparticles and hence the lifetime and branching ratios for possible decay modes of the 
gravitino depend on the mass spectrum. We calculate these quantities in detail, as we 
explain below. 

In order to calculate the decay rate of the gravitino, it is necessary to fix the mass spec- 
trum and the mixing matrices in the MSSM sector. Although the effects of the gravitino 
on the BBN can be calculated for arbitrary mass spectrum of the MSSM particles, it is not 
practical to study all the possible cases since there are very large number of parameters 
in the MSSM sector. Thus, we adopt a simple parameterization of the SUSY breaking 
parameters, that is, the mSUGRA-type parameterization of the soft SUSY breaking pa- 
rameters. We (parameterize the MSSM parameters by using unified gaugino mass mi/2, 
universal scalar mass mo, universal coefficient for the tri- linear scalar coupling Aq, ratio of 
the vacuum expectation values of two Higgs bosons tan (3, and supersymmetric Higgs mass 
l^H- Then, the properties of the superparticles (including the gravitino) are determined 
once these parameters as well as the gravitino mass ^3/2 are fixed and, consequently, 
we can derive the upper bound on the reheating temperature. Notice that, although we 
adopt the simple parameterization of the soft SUSY breaking parameters, our analysis is 
applicable to more general cases as far as the gravitino is heavier than one of the MSSM 
superparticles (like the lightest neutralino). 

Even with the mSUGRA parameterization of the MSSM parameters, the whole pa- 
rameter space is still too large to be completely studied. Thus, in this paper, we pick up 
several typical mSUGRA points and derive the constraints for these points. In particular. 
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Case 1 Case 2 Case 3 Case 4 

mi/2 300 GeV 600 GeV 300 GeV 1200 GeV 

mo 141 GeV 218 GeV 2397 GeV 800 GeV 

Aq 

tan/3 30 30 30 45 

liH 389 GeV 726 GeV 231 GeV -1315 GeV 

m^o 117 GeV 244 GeV 116 GeV 509 GeV 

^(thcrmal)^2 q^^^^ q^^^q q^^qq Q.lll 

Table 2: mSUGRA parameters used in our analysis. 



we pick up points where the thermal rehc density of the LSP becomes consistent with the 
dark matter density determined by the WMAP observation The points we consider 
are hsted in Table El For all the cases, we checked that the lightest neutralino becomes 
the LSP (if the gravitino mass is larger than Tn^o). Using the DarkSUSY package [23], we 

calculated the thermal relic density of the LSP Q^^p^'^^\*^ (We use h = 0.71 jl], where h 
is the expansion rate of the universe in units of 100 km/sec/Mpc.) 



3 Gravitino Decay 

In order to study the effects of the gravitino on the BBN, it is important to understand 
the lifetime of the gravitino and its decay modes. Since the gravitino is the gauge field 
for the supersymmetry, gravitino couples to the supercurrent and hence its interaction is 
unambiguously determined. Possible decay modes have, however, model-dependence. In 
the following, we discuss how the decay rate and the branching ratios of the gravitino are 
calculated. 



3.1 Interaction and two-body processes 

We first briefly summarize the interactions of the gravitino.*^ Gravitino is the super- 
partner of the graviton, and it couples to the supercurrent. Thus, (relevant part of) the 
interaction of the gravitino is given by 



G 



E [x^^'VYPLi'.D^^P^^^ + h.c] , (3.1) 



V2M, c 



#5The Cases 1 and 2 are in the so-called "co-annihilation region," while the Cases 3 and 4 are in the 
"focus-point region" and "Higgs funnel region," respectively. 
T^^For details, see, for example, 
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(G) 



where the sum over G is for all the gauge multiplets (consisting of the gauge field A 
and the gaugino A*-*^^) while C for the chiral multiplets (consisting of the scalar boson 0*^*^^ 



and the fermion x'^'^^)- Here, -Fp^-* is the field strength for A^^\ and D^, represents the 
covariant derivative. In addition, ~ 2.4 x 10^^ GeV is the reduced Planck scale. As 
is obvious from Eq. (jSUl), if we restrict ourselves to consider the two-body final states, 
gravitino decays into some standard-model particle and its superpartner. 

From the Lagrangian given in Eq. (j3.1|) . we can read off the vertex factors for the 
gravitino and calculate the (partial) decay rates of the gravitino. We first consider the 
decay processes with two-body final states. Then, the decay rate is expressed as 



(G) 



/3fiVc 



gauge 



167rm3/2M2 



X \M\\ 



(3.2) 



where M. represents the Feynman amplitude for the decay process (with = 1), and A^^c 
is the color factor: Nq = 8 for the process ipfi 99, = 3 for the processes with quark 
and squark final state, and Nq = 1 otherwise. In addition, for the process '^p^i ^B, /?f 
is given by 



1 



2 N211/2 



m 



3/2 



- 2m2/2(m^ + m|) + (m^ - m|)2] 



(3.3) 



where and are masses of A and B, respectively. 

For the gravitino decay process into a gauge boson V and a fermion x, we define p, q, 
and q' to be the momenta of tp^, V, and respectively. Then, we obtain 



pq 



pq 



qq 



("^3/2 + - ml) 

{ml/2 -rnl + mj) 
mz) 



2 ("^3/2 



m 



V 



2' 

x) ' 



(3.4) 
(3.5) 
(3.6) 



where my and are masses of V and Xi respectively. With these quantities, for the 
process with a massless gauge field in the final state (i.e., iIj^ —>■ 7%° and gg), we obtain 



\M\' 



1,9 



- cr'c 



(G)^(G)* 



{pqf{pq') 



m 



3/2 



(3.7) 
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while with massive gauge field in the final state (i.e., ipfj, Zx^i and W^xf) > we obtain 
2 



\M\' 



MG)MGh , ^(G)^(G) 
(G)^(G)* 



2 + m)m ) - TS'^viPQ ) 



m 



3/2 



(^)^W* + h.c.) ms/2 
+ (cfcr-Cfcr + h.c.)m,(,g) 



m 



3/2 



.(H)^(H)h 



+3 (crc^ 



(H)^(H)* 



^(H)^(H)h 



1 + 



2"^3/2"^V' 



Here, Cf \ and C)f^ are vertex factors. These vertex factors for individual 

processes are given in Appendix El For the decay processes with a scalar boson in the final 
state, we identify p, g, and g' to be the momenta of t/'^, scalar boson 0, and fermion x, 
respectively. Then, the products of the momenta can be obtained from Eqs. ()3.4|) — ()3.6|) 
by replacing my — > m^^ with being the mass of the scalar boson. Then, we obtain 



^(H) 



^(H) 



2 ^ 1 

scalar o 



{pqf 



m 



3/2 



^(C)^(C)* 
L L 



The vertex factors C^^^ and Cj^"* are also given in Appendix El 



(3.9) 



.(C) 



3.2 Three-body processes 

In this paper, we consider the case where the LSP is the lightest neutralino Xv Then, the 
two-body process ipfj, — >• 7X1 is always allowed, and the (total) decay rate of the gravitino 
is determined by the two-body process(es). In studying the effects of the gravitino on the 
BBN, however, it is also important to understand the spectrum of the hadrons produced 
by the decay of the gravitino. 

In most of the cases, the number of the hadrons from the gravitino decay is mostly 
determined by the two-body processes. For example, if the gravitino can decay into 
a superparticle other than the LSP, the emitted superparticle decays into Xv this 
secondary decay process, sizable number of the hadrons may be produced. In addition. 
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when the mass difference between the gravitino and the hghtest neutrahno is larger than 
mz, the decay process ~^ becomes kinematically allowed. In this case, the decay 
of the Z-boson produces large amount of the hadrons. In most of the cases, the number 
of the hadrons produced from those two-body processes is much larger than that from the 
three body processes. Then, the three-body processes are irrelevant for our study. 

However, in some case, precise determination of the hadron spectrum requires the 
calculation of the three-body final state processes. In particular, the three-body processes 
become important if (i) 777.3/2 — m^o < mz, and (ii) all the superparticles except Xi and 
sleptons are heavier than the gravitino. Notice that, when the condition (i) is satisfied, 
it may be the case that the only possible two-body decay process is ipfj, ■jXi hence 
the hadrons are not produced by the two-body process. In some case, gravitino may 
also decay into lepton and slepton pair, but the decays of the (s)leptons does not produce 
significant amount of hadrons in our case. Thus, we pay particular attention to the process 
qqXi when 7773/2 — 777^0 < 777^. 

When the LSP is the neutralino, the relevant three-body processes are induced by the 
diagrams shown in Fig. ^ In our study, we consider the effects of all the diagrams listed 
in Fig. [Hand calculate the decay rate for the process — > qqXi- To see the importance 
of the 3-body process, in Fig. |21 we plot the "3-body hadronic decay width" defined as*^ 

^{^, ^ qqxl) = ^{'4^, ^ uuxl) + r(^^ ^ ddx\) + r(^^ ^ ssxl) 

+r(7/;^ ^ ccxl) + r(7/;^ ^ hhxl) + r(7/;^ ^ ttx\). (3.10) 

In Fig. the MSSM parameters are taken to be those for the Case 1. In this case, the 
3-body decay is induced dominantly by the photon-mediated diagram. When the mass 
difference between the gravitino and the LSP becomes larger than 777^, however, Z-boson 
mediated contribution with ~ mz becomes the most important, where m^q is the 
invariant mass of the qq system. In fact, such a process should rather be classified into 
the 2-body process ipfj, — > Zxi followed hy Z qq. To see this more explicitly, we also 
plot the quantity r(7/^^ Zxi) x Br(Z — > qq). As one can see, T{ip^j, — ^ qqXi) is well 
approximated by r(?/^^ — >• Zxi) x Br(Z — >• qq) when the decay process V'/i Zxi becomes 
kinematically allowed. 

In our numerical study, we treat the process t/)^ qqXi the following way: 

• When 7773/2 — 777^0 < mz, we calculate the Feynman diagrams shown in Fig. ^and 
calculate the decay rate r(7/;^ qqXi)- (J^ this case, the decay process 7/;^ — > Zxi 
is kinematically forbidden and hence is irrelevant.) 

• When 7773/2 — 777^0 > 777^, we approximate the hadronic decay rate induced by the 
diagrams in Fig. ^by T^ipf^ — * Zxi) x Bt{Z qq). 

we discuss in the following, the three-body process becomes important when the mass difference 
between the gravitino and the LSP is smaller than mz- Thus, the process 1/'^ — > ttxi ^'^^ important for 
our analysis. 
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(c) 



(d) 



Figure 1: Feynman diagrams for the process ip/j 



The "blobs" are from the 



gravitino-supercurrent interaction. For (d), there is also CP-conjugated diagram (with 
the replacements q ^ q and q* q)- 



With our treatment, the effect of the photon-mediated diagram (as well as those from 
Figs, life) and (d)) is neglected when m^i2 — m^o > rriz- However, effect of such diagram 
is subdominant since the process is mainly induced by the Z-boson mediated diagram. 

When 1713/2 — m^o > mz, energy distribution of the quark and anti-quark is easily 
calculated since the decay is dominated by the process with nigg ~ mz- When 1713/2— m^o < 
mz, on the contrary, mqq has broader distribution. Thus, in this case, we numerically 
calculate the differential decay rate 

dTjipf, qqxi) 

to obtain the energy distributions of the quarks and anti-quarks emitted from the three- 
body processes. (In the calculation of dV^ipfj, —>■ qqXi) / dm'^g, we approximated that the 
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20 40 60 80 100 120 



m3/2 - m^o (GeV) 

Figure 2: Width for the process — >• qqXi as a function of 1713/2 — m^o (sohd hne). 
We adopt the mSUGRA parameters for the Case 1. For comparison, the decay rate 
r('?/'/^ —>■ Zxi) X Bt{Z —>■ qq) is also shown in the dashed hne. 



final-state q and q have isotropic distribution in their center-of-mass frame.) The hadron 
spectrum from the three-body decay process is obtained by using this differential decay 
rate. When 171^/2 — m^o < mz, the photon-mediated diagram gives the dominant contri- 
bution while the effects of other diagrams are almost irrelevant (unless 712^/2 — m^o is very 
close to mz). In Appendix El we present the approximated formula of the differential 
decay rate, only taking account of the photon-mediated diagram. 

3.3 Lifetime and branching ratios of the gravitino 

Now we can quantitatively discuss the decay rates of the gravitino. First, following the 
procedures discussed in the previous subsections, we calculate the partial decay rates of 
the gravitino for all the possible decay modes. Adding all the contributions, we obtain the 
lifetime of the gravitino: 

^3/2 = T^T—^ TfT- (3.11) 

T[i)^^ all) 

We calculate r3/2 as a function of the gravitino mass for the cases listed in Table El and 
the results are shown in Fig. El Importantly, lifetime of the gravitino becomes shorter as 
the gravitino becomes heavier. 

As one can see, when the gravitino mass is smaller than a few TeV, r3/2 for the Case 4 
is found to be longer than those for other cases. This is due to the fact that, for the Case 
4, masses of the MSSM particles are larger than other cases and hence the decay rates 
of the gravitino in this case is suppressed. When the gravitino is much heavier than the 
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MSSM particles, on the contrary, the hfetime of the gravitino is insensitive to the mass 
spectrum of the MSSM particles. 

Importantly, when the gravitino is lighter than ~ 20 TeV, r3/2 becomes longer than 
1 sec and hence the relic gravitinos decay after the BBN starts. Thus, in particular in this 
case, significant constraints on the reheating temperature after the inflation is expected. 

Branching ratios of the decay process of the gravitino depends on the model parameters. 
To see this, for the Cases 1 — 4, we plot the branching ratios for various two-body final 
states in Fig. |3 As one can see, the branching ratios have sizable model dependence when 
the gravitino mass is relatively small. This is because, when the gravitino mass is small, 
decay rate of the gravitino is sensitive to the mass spectrum of the superparticles. When 
m3/2 becomes large enough, on the contrary, branching ratios become insensitive to the 
model parameters. 

4 Secondary Decays and Hadronization 

Although the gravitino primarily decays into a standard-model particle and its superpart- 
ner (or into the qqXi final state), most of the daughter particles also decay with time scale 
much shorter than the cosmological time scale. In addition, all the partons (i.e., quarks 
and gluon) are hadronized into mesons or baryons. These processes are important in the 
study of the BBN with the primordial gravitinos and, in this section, we discuss these 
issues. 

The possible decay modes of the individual superparticles strongly depend on the mass 
spectrum as well as on the mixing and coupling parameters. In order to systematically 
take account of all the relevant decay processes, we use the ISAJET package jSSj which 
automatically calculates the partial decay rates of all the unstable particles. 
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Figure 4: Branching ratios of the decay of the gravitino as functions of the gravitino mass. 
The thick sohd hne is for the final states Xi + anything, dot-dashed fine for lepton-slepton 
pairs, dotted hne for Xi = 2 — 4) or chargino + anything, dashed hne for gluon-gluino 
pair, and thin sohd hne for quark-squark pair final states. 



In order to discuss the hadro-dissociations induced by the gravitino decay, we should 
first calculate the spectra of the partons (i.e., u, d, s, c, b, t and their anti-particles, and 
gluon). There are two types of processes producing energetic partons: one is the decay of 
the gravitino and the other is the subsequent decays of the daughter particles. Spectra of 
the partons of the first type are directly calculated by using the partial decay rates of the 
gravitino presented in the previous section. (Here, the effects of the "3-body" processes 
are also included when 7713/2 — rn^o < mz-) In order to calculate the contribution of the 
second type, we have to follow the decay chain of the unstable particles. In addition, 
the parton spectra should be calculated by averaging over all the possible decay processes 
with the relevant branching ratios of individual particles. In our analysis, the decay chain 
is followed by using the PYTHIA package [20! which automatically take account of the 
decay processes of the unstable particles (including the superparticles). 

At the cosmic time relevant for the BBN, time scale for the hadronization is much 
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Ekm (GeV) Ekm (GeV) 

Figure 5: Distributions of the nucleons (i.e., (a) proton and (b) neutron) from the decay of 
a single gravitino as functions of the kinetic energy. Here we take the mSUGRA parameters 
for the Case 1, and m3/2 = ITeV. 



shorter than the time scale for the scattering processes and hence all the partons are 
hadronized before scattering off the background particles. Thus, it is necessary to calculate 
the spectrum of the mesons and baryons produced from the partons. In particular, for our 
analysis of the BBN-related processes, proton, neutron, and charged pions play significant 
roles. 

In our analysis, the hadronization processes are dealt with the PYTHIA package fl^ : 
we have modified the PYTHIA package to include the primary decay processes of the 
gravitino, then the subsequent decay processes of the daughter particles (including the 
superparticles) and the hadronization processes of the emitted partons are automatically 
followed by the original PYTHIA algorithm. In Fig. El we plot the distribution functions 
of the proton and neutron as functions of their kinetic energy (i.e., i?kin = E — 
with m^v being the corresponding nucleon mass). In order to check the reliability of our 
estimate of the hadron spectrum, we have performed an independent calculation using the 
ISA JET package [25j; we have also modified the ISA JET package to include the decay 
processes of the gravitino and we calculated the hadron spectrum. We have checked that 
the difference between the two calculations is within 10 % level. In particular, for the 
region with relatively large kinetic energy, which gives the most important effects on the 
hadro-dissociation processes of the light elements, difference is very small. 

Before closing this section, we define one important parameter, which is the (averaged) 
visible energy emitted from the gravitino decay. Once a high energy particle with electro- 
magnetic interaction is injected into the thermal bath, it induces the electro-magnetic 
shower and, consequently, the photo-dissociation processes of the light elements are in- 
duced by the energetic photons in the shower. The event rates of the photo-dissociation 
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processes (for a fixed background temperature) are mostly determined by the total amount 
of the "visible" energy injected into the thermal bath jB]. Since we consider the case where 
the i?-parity is conserved, some fraction of the energy is always carried away by the LSP 
for the decay process of the gravitino. Taking account of such effect, we calculate the 
averaged visible energy emitted by a single decay of the gravitino:*® 

^vis = ^3/2 - (^^o) - (E,), (4.1) 

where the second and third terms of the right-hand side are the (averaged) energy carried 
away by the LSP and the neutrinos, respectively. -Evis is used for the calculations of the 
photo-dissociation rates. 



5 Light-Element Abundances 
5.1 Theoretical calculation 

Now, we explain how we calculate the light-element abundances. In order to set a bound on 
the reheating temperature after the inflation, we assume that the gravitinos are produced 
by the scattering processes of the thermal particles. Using the thermally averaged gravitino 
production cross section given in [211, the "yield variable" of the gravitino, which is defined 
as F3/2 = is given by jTE] 



3/2 



1.9 X 10 



-12 



R 



IQio GeV 



1 + 0.045 In 



1010 GeV 



1 - 0.028 In 



R 



1010 GeV 



(5.1) 



where 713/2 is the number density of the gravitino while s = ^g^,s{T)T^ is the entropy 
density with (7^,5 (T) being the effective number of the massless degrees of freedom at the 
temperature T, and the reheating temperature is defined as*^ 



R 



10 



2 * int 



1/4 



(5.2) 



with Finf being the decay rate of the infiaton. 

'^^For unstable leptons and mesons (in particular, pions), we checked that their lifetimes are shorter 
than their mean free time. Thus, in the calculation of the visible energy, we treated them as unstable 
particles and hence the energy carries away by the neutrinos are not included in i?vis- 

'^^Strictly speaking, the "reheating temperature" corresponds to the maximal temperature of the last 
radiation dominated era. Thus, if some scalar field (f) other than the inflaton once dominates the universe 
after the inflation, the reheating temperature here is given by the same expression as Eq. (|5.2|) with Finf 
being replaced by the decay rate of <j>. One of the examples of such scalar fields is the curvaton |28| which 
provides a new origin of the cosmic density perturbations. 
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Once produced, the gravitinos decay with very long hfetime. In particular, if the 
gravitino mass is smaller than ~ 20 TeV, gravitinos decay after the BBN starts and hence 
the light-element abundances may be significantly affected. In order to study the BBN 
processes with unstable gravitino, our study proceeds as follows: 

1. MSSM parameters are determined for one of the mSUGRA points given in Table |21 
Then, all the mass eigenvalues and mixing parameters are calculated. 

2. Using the parameters given above, we calculate partial decay rates of the gravitino 
for all the kinematically allowed 2-body processes. When 7713/2 —^xi ^ 
calculate T{ip^ QQXi)- 

3. We perform a Monte-Carlo analysis using the branching ratios of the gravitino to 
calculate the energy distribution of the hadrons produced by the decay of the grav- 
itino. As explained in the previous sections, the decay chain of the decay products 
(including the superparticles) and the hadronizations of the emitted partons are fol- 
lowed by the modified PYTHIA code. Simultaneously, we calculate the (averaged) 
emitted visible energy from the decay of the gravitino. 

4. For a given reheating temperature Tr, we calculate the abundance of the thermally 
produced gravitino using Eq. (jS.lj) . 

5. We calculate the light-element abundances, taking account of the decay of the ther- 
mally produced gravitinos. (Standard reactions of the light elements are also in- 
cluded.) We use the baryon-to-photon ratio suggested by the WMAP ji]: 

7/ = (6.1 ±0.3) X 10"^°, (5.3) 

at the la level. (Here we enlarged the lower error bar from 0.2 to 0.3 since the 
Monte-Carlo simulation is easier if the error bar is symmetric. This does not signif- 
icantly change the resultant constraints.) The baryon-to-photon ratio is related to 
the density parameter of the baryon as flsh'^ = 3.67 x lO^r/. 

6. Since the event rates of the non-standard processes induced by the gravitino decay 
are proportional to the abundance of the primordial gravitinos, deviations of the light 
element abundances from the standard BBN results become larger as the reheating 
temperature becomes higher. The resultant light element abundances are compared 
with the observational constraints and upper bound on the reheating temperature 
is obtained. 

Although the details of the analysis of the photo- and hadro-dissociation processes 
and p ^ n interchange are explained in El El ^] , we briefly summarize several im- 
portant points. Once energetic hadrons are emitted into the thermal bath in the early 
universe, they induce the hadronic shower and energetic particles in the shower cause 
hadro-dissociation processes. In addition, released visible energy from the gravitino decay 
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Figure 6: Number of hadrons(7r, p, n) produced by the decay of single gravitino for the 
Case 1. (The number of vr" is the same as that of vr^ within the expected error of the 
Monte-Carlo analysis.) 



eventually goes into the form of radiation which cause electro-magnetic shower. Energetic 
photons in the shower cause photo- dissociation processes. Furthermore, when the cosmic 
temperature is relatively high (T^ 0.1 MeV), p <-> n inter-converting processes by the 
nucleons and the charged pions become significant. 

When the cosmic temperature is higher than 0.3 MeV, the p ^ n inter-converting 
processes induced by the charged pions (i.e., p + vr" ^ n + and n + tt~^ ^ p + 7r°) are 
the most important. Since the resultant ^He abundance is sensitive to the n/p ratio, such 
inter-converting processes affects the ''He abundance. 

Since the charged pions are expected to be stopped in the thermal bath before inter- 
converting the background nucleons, we need to know only the total numbers of 7r+ and 
TT^ produced by the decay of the gravitino. The number of pions produced by the decay 
of the single gravitino is plotted in Fig. IHl as a function of the gravitino mass for the Case 
1. As one can see, the number of the pions increases as the gravitino mass becomes larger. 
In addition, when the gravitino mass is smaller than ~ 1 TeV, partial decay rates of the 
gravitino have significant dependence on m3/2 because the gravitino mass becomes close 
to the masses of the MSSM superparticles in this region; consequently, number of pion 
has strong dependence on 171^/2. In the figure, we also plotted the number of the proton 
and the neutron produced by the decay of the gravitino. 

As the lifetime of the gravitino becomes longer, it is likely that most of the thermally 
produced gravitinos decay after ^He and other light elements (like D, T, '^He, and so on) are 
synthesized. Then, the hadro- and photo-dissociation processes may significantly change 
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the light-element abundances. 

When 10^ sec ^ r3/2 ^ 10^ sec, hadro- dissociation processes of the light elements are 
important while, for longer lifetime, photo-dissociation processes give more significant 
constraints. In particular, since the number density ^He is much larger than those of other 
light elements, dissociation of ^He may significantly change the abundances of D and '^He. 
In addition, non-thermally produced T and ^He may scatter off the background ^He to 
produce ^Li via the processes T + ^He ^Li + p and ^He -|- ^He — > ^Li + n jTUl CH 113 Hn] . 
Since there is very stringent upper bound on the primordial abundance of ^Li, such non- 
thermal production process of ^Li gives significant constraint on the reheating temperature. 

5.2 Observational constraints 

In order to derive a bound on the reheating temperature, we compare the theoretical results 
of the light-element abundances with the observational constraints. The observational 
constraints we use are summarized below. Since there are uncertainties in the constraints, 
we adopt several different bounds on the primordial abundances of the light elements 
in some case. The errors of the following observational values are at la level unless 
otherwise stated. When we adopt both the statistical and systematic errors, we add them 
in quadrature. 

For the D abundance, we use constraints obtained from the measurements in the high 
redshift QSO absorption systems [SH]. Here, we consider two constraints; one is the 
"averaged" constraint 

K/nH)""^ = (2.78+°:^^) X 10-^ (5.4) 

while the other is the "conservative" one, which is the highest value of UD/nu among the 
results listed in 

(no/nH)""^ = (3.98+°:f7) x 10"^ (5.5) 

(Here and hereafter, the superscript "obs" is used for the primordial values inferred by 
the observations.) 

Abundance of ^He may significantly change during the evolution of the universe from 
the BBN epoch to the present epoch. Thus, for ^He, it is not easy to observationally 
determine its primordial value. In our analysis, we do not rely on any detailed model of 
chemical evolution to derive bound on the primordial abundance of ^He. Instead, we only 
use the fact that D is more fragile than '^He. Then, we expect that the ratio rs 2 = n3He/''^D 
does not decrease with time jHOl [13 CHI EI- The solar-system value of ^He-to-D ratio is 
measured as 

= 0.59 ±0.54 {2a). (5.6) 
Thus, we obtain the upper bound on the primordial ^He to D ratio 

rt2 < rf,2- (5.7) 

For the primordial abundance of ^He, we use the constraints from the recombination 
lines from the low metallicity extragalactic HII regions. Taking into account the fact 
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that several groups independently derived bounds on the mass fraction of "^He, we derive 
upper bound on Tr with the following three different bounds: the first one is based on the 
analysis by Fields and Olive jHSj 



F°bs(pQ) _ Q 238 ± (0.002)stat ± (0.005)syst, (5.8) 

where the first and second errors are the statistical and systematic ones, respectively, the 
second is obtained by Izotov and Thuan |38] : 

F°^^(IT) = 0.242 ± (0.002),tat(±(0.005)syst), (5.9) 

where we have added the systematic errors following Refs. [311E3ISni, and the last one is 
by Olive and Skillman [87j : 

r°bs(Qg^) _ Q ± 0.009, (5.10) 

where the error includes both the statistical and systematic ones. With these three con- 
straints, we will discuss how the upper bound on Tr changes as we adopt different value 
of 

The primordial value of the ''Li abundance is observed in old Pop II halo stars. Typi- 
cally riTLi/riH is O(10~^°). In [31] , it was reported that 

(mLi/nH)°'^ = 1.23l°J^ X 10-10, (5.11) 

while, recently, relatively higher value of the ''Li abundance was also reported ^U^: 

log.olinru/nnT'''] = -9.66 ± (0.056)stat ± (0.06)syst. (5.12) 

Here, '''Li abundances given in 1)5.111) and ()5.12)1 differ by the factor ~ 2, and we expect 
that there is still some large uncertainty for the observational values of mu/nu- Here, we 
adopt the constraint given in [3^] with an additional large systematic error, considering 
the possibilities of the increase by the cosmic-ray spallation of the C, N, O and so on, and 
the decrease by the depletion by the convection in the stars jlO] 



log 



10 



(^Li/nn)"^' = -9.66 ± (0.056),tat ± (0.300)add. (5.13) 



The linear-scale value is given by {nTu/nn)"^^ = (0.54 — 8.92) x 10^^° at the 2a level. In 
deriving the upper bound on Tr, we use the constraint (|5.13|) . 

Usually the abundance of ^Li is measured as a ratio of ^Li and ^Li in the old Pop II halo 
stars; at the 2a level, {nau/mu)^'^^" = 0.05 ± 0.02 [H]. The primordial value is expected 
to be smaller than this value because it is likely that the the cosmic-ray spallation has 
produced additional ^Li after the BBN 1 ^ Elj. By adopting the milder value of the 
constraint on mu/nu given in Eq. ()5.13|) . we get the upper bound on the primordial value 
of nen/nn at the 2a level, 

(n6Li/nH)°"^ < (l-lOtl-lt) X 10-11 (2a). (5.14) 
We use this value as the upper bound on neu/na except in Section IFT^ 
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6 Numerical Results 



6.1 Upper bound on Tr 

Now, we are at the position to show our numerical results. As discussed in the previous sec- 
tions, we calculate the light element abundances as functions of the gravitino mass, other 
MSSM parameters, and the reheating temperature Tr. Then, we compare the theoretical 
prediction with the observations. In order to systematically derive the upper bound, we 
calculate the variable defined as 

= (^th)2 ^ ('^obs)2 for X, = (nD/nH),y, (6.1) 

where x^^ and x°^^ are the center values of Xi determined from the theoretical calculation 
and observations, while af^ and af^^ are their errors, respectively. In our analysis, (ct|^)^ 
is calculated by the Monte-Carlo analysis. For Xi = 2 {nei^-Jnu) and log]^Q[(?27Li/nH)] we 
only use the upper bound, and we define Xi as*^° 

{('™th ™obs\2 
) fora;, = r3,2, (n6Li/nH),logio[(mLi/nH)]. (6.2) 

: otherwise 

With these quantities, we derive 95 % level constraints which correspond to xf = 3.84 for 
Xi = (no/nn) and Y, and xl = 2.71 for Xi = r3,2, {neu/rin) and logioKn? u/nn)]- (For 
details, see [IE]-) 

In Figs. |71 — Uni we show the upper bounds on the reheating temperature. For D, 
we considered the observational constraints ()5.4p and ()5.5|) to see how the upper bound 
depends on the bound on D. For ^He, ''Li, and ^Li, we use fl5.6|) . ()5.13|) . and ()5.14|) . 
respectively. For ^He, we consider three cases fjS.Sj) (FO), ()5.9p (IT), ()5.10p (OS), since the 
upper bound on Tr for the case of relatively heavy gravitino is sensitive to the observational 
constraint on the abundance of ^He. In deriving Figs. [7| — ^1 the MSSM parameters are 
determined by using the mSUGRA parameters given in TableEl In addition, the lifetime of 
the gravitino as well as its branching ratios are calculated using the MSSM mass spectrum 
obtained from these parameters. In this analysis, we concentrated on the case where the 
gravitino is unstable. In the figures, we shaded the region where m3/2 < "^xl- 

Although we have considered four different cases with different mass spectrum of the 
MSSM particles, the qualitative behavior of the constraints are quite insensitive to the 
choice of underlying parameters. When the gravitino mass is larger than a few TeV, most 
of the primordial gravitinos decay at very early stage of the BBN. In this case, in addition, 

#iOpQj- riTLi, we only use the upper bound since we could not include one of the non-thermal production 
process of mi^^, because of the lack of the experimental data. For details, see and also the next 
subsection. 
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Figure 7: Upper bound on the reheating temperature for the Case 1 as a function of the 
gravitino mass. 



photo- and hadro-dissociations are ineffective. Then, overproduction of ^He due to the 
p n conversion becomes the most important. For the observational constraints on the 
mass fraction of ^He, we consider three different observational results given in Eqs. fj5.8|) — 
f)5.10|) . As one can see, the upper bound on Tr in this case is sensitive to the observational 
constraint on the primordial abundance of ^He; for the case of ^3/2 = 10 TeV, for example, 
Tr is required to be lower than 3 x 10^ GeV if we use the lowest value of Y given in Eq. 
(I5.8|l while, with the highest value given in Eq. ()5.1()j) . the upper bound on the reheating 
temperature becomes as large as 4 x 10^ GeV. 

When 400 GeV ^ 777.3/2 ^ 5 TeV, gravitinos decay when the cosmic temperature is 1 keV 
— 100 keV. In this case, hadro-dissociation gives the most stringent constraints; in particu- 
lar, the overproductions of D and ^Li become important. Furthermore, when the gravitino 
mass is relatively light (m3/2^400 GeV), the most stringent constraint is from the ratio 
■^He/D which may be significantly changed by the photo-dissociation of ^He. 

It should be noted that, even when the gravitino cannot directly decay into colored 
particles (i.e., the squarks, gluino, and their superpartners) due to the kinematical reason, 
the reheating temperature may still be stringently constrained from the hadro-dissociation 
processes. This is due to the fact that some of the non-colored decay products (in partic- 
ular, the weak bosons and Z as well as some of the superparticles) produce hadrons 
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Figure 8: Same as Fig. [7| except for the MSSM parameters are evaluated for the Case 2. 



when they decay. In particular, when the mass difference m^i2 — m^o is larger than m^, 
hadro-dissociations become important since sizable amount of hadrons are produced by 
the process ipfj, ~^ ^Xi- 

Although Figs. [7| — [TUl look roughly the same, the upper bound on the reheating 
temperature has model dependences. In particular, for a fixed value of the gravitino mass, 
the upper bound depends on the MSSM parameters as one can see in the figures. To see 
this in more detail, in Table El we show the upper bound on Tr for the cases listed in 
Table El For the fixed value of the gravitino mass, the upper bound on Tr may vary by 
the factor as large as ~ 10 when the gravitino mass is of the order of 1 TeV. When the 
gravitino becomes heavier than ~ 10 TeV, however, the upper bound becomes insensitive 
to the model parameters. This is due to the fact that, in this case, the branching ratios 
of the gravitino are almost independent of the MSSM parameters. 

Model-dependence of the upper bound on the reheating temperature is mostly from 
the change of the lifetime and decay modes. For the Case 1, we have chosen the mSUGRA 
parameters which give relatively light superparticles (compared to other cases). In the 
Case 2, masses of all the superparticles are slightly increased compared to the Case 1. 
Consequently, we can see some changes of the constraints on the reheating temperature. 

Compared to the Case 1, scalar masses are significantly increased in the Case 3 with 
the gaugino masses being unchanged. In this case, gravitino is likely to decay into the 
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gauginos, in particular, into the gluino when kinematically allowed. (See Fig. |3) We 
found that the gluon-gluino final state produces more hadrons (in particular, protons and 
neutrons) than the quark-squark final state. Consequently, in the Case 3, upper bound on 
Tr becomes lower than that for the Case 2. We have also studied the case where masses 
of all the squarks and sfermions are pushed to infinity by hand while keeping the gaugino 
mass as low as 0(100 GeV). In this case, the constraint on Tr is almost the same as that 
for the Case 3. In addition, in the Case 4, masses of all the superparticles are very large 
(~ a few TeV). Then, lifetime of the gravitino becomes relatively long, which makes the 
upper bound less stringent for gravitinos with 777.3/2 ~ a few TeV. 

Although our main concern is to study the effects of the gravitino decay on the BBN, 
it is also important to consider other constraints. One of the important constraints is 
from the production of the LSP from the decay of the gravitino. Importantly, the LSP 
is produced with the decay of the gravitino, and the present number density of the LSP 
is given by the sum of two contributions; thermal relic, which is calculated with the 
DarkSUSY package for each cases, and the non-thermally produced LSP from the gravitino 
decay. Since one LSP is produced by the decay of one gravitino, the density parameter of 
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the LSP which has non-thermal origin is given by*^^ 

An,,pft==. 0.054 x(^)(^^). (6,3) 

where we have neglected the logarithmic corrections in Eq. (jS.lj) . If we require, for example, 
that the total mass density of the LSP be within the 95 % CL. bound of the WMAP 
constraint (i.e., ^l^^p^^^^h"^ + Af^LSP^^ < 0.1287 we also obtain upper bound on Tr, 
which is given by 3 x 10^ GeV, 1 x 10^ GeV, 4 x 10^ GeV, and 6 x 10* GeV, for the Cases 
1, 2, 3, and 4, respectively.*^^ In our numerical analysis, we calculated the abundance of 
the LSP taking account of the entropy production by the decay of the gravitino; constraint 
from Af^LSP is also shown in the figures. 

Another constraint may be obtained from the distortion of the cosmic microwave back- 
ground (CMB). An additional injection of the photon into the thermal bath by the decaying 
particles is severely constrained in order not to disturb the black-body shape of the CMB 

'^^^In fact, entropy production occurs when the gravitino decays, and consequently, primordial LSPs are 
diluted by some amount. For the reheating temperature giving the constraint on the relic density of the 
LSP, however, the effect of the dilution is negligible. 

'^^^This bound is sensitive to the choice of the MSSM parameters since the abundance of the thermally 
produced LSPs depends on the MSSM parameter. If we choose a parameter set which gives fi^^gp™'*''' 
much smaller than the WMAP value, bound from the production of the LSP may become much weaker. 
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Case 1 



Case 2 



Case 3 



Case 4 



300 GeV 

1 TeV 

3 TeV 
10 TeV (FO) 
10 TeV (IT) 
10 TeV (OS) 



6 X 10^ GeV 
5 X 10^ GeV 
1 X 10^ GeV 

3 X 10^ GeV 
8 X 10* GeV 

4 X 10^ GeV 



3 X 10^ GeV 
1 X lO*' GeV 
1 X 10^ GeV 

3 X 10^ GeV 
8 X 10* GeV 

4 X 10^ GeV 



4 X 10^ GeV 
4 X 10^ GeV 
4 X 10^ GeV 

2 X 10^ GeV 
6 X 10* GeV 

3 X 10^ GeV 



6 X 10^ GeV 
1 X 10^ GeV 

3 X 10^ GeV 
8 X 10* GeV 

4 X 10^ GeV 



Table 3: Upper bound on Tr for several values of the gravitino mass. For the D, '^He, and 
^Li abundances, here we use the observational constraints ()5.4|) . ()5.6|) . ()5.14p . For ^He, we 
consider three cases dHU) (FO), dHSD (IT), (UnUI) (OS). 



spectrum |45| I46j : numerically, \fi\ < 9 x 10^^ for yU-distortion and \y\ < 1.5 x 10~^ for 
y-distortion are required [Hj. Using these constraints, we obtain the upper bound on the 
total amount of the injected energy Ap^; using the relation Ap^/s = £^vis^3/2, 

^ < 1.60 X 10--GeV X (^) ''^^ exp [{r,c/ry,f'] , (6.4) 

for /i-distortion for rjc ^ T3/2 ^ 2.5 x 10^ sec x (fife/i^/0.022) with Tdc being decoupling 
time of the double Compton scattering: 

Tdc = 6.10 X lO^sec X — — ^ —] , 6.5 

^ V2-725 K J VO-022 J \ 0.88 / ' ^ ' 

and 

^<2.7xlO-3GeVxf^^)"'', (6.6) 
s V 1010 sec/ ' ^ ^ 

for ^/-distortion for Ts/2 ^ 2.5 x 10^ sec x (fi;,/i^/0.022) jl^. Here To is the photon temper- 
ature at present. Constraint from the distortion of the CMB spectrum is also shown in 
the figures. 



6.2 Comment on the Li abundance 

So far, we have considered the constraints on the reheating temperature, assuming that 
the prediction of the standard BBN agrees with observations. Although the standard 
BEN predicts the light-element abundances which are more or less consistent with the 
observational constraints, however, it has been pointed out that, if we adopt the baryon-to- 
photon ratio suggested by the WMAP, standard BBN predicts the ^Li abundance slightly 
larger than the observed value. Indeed, if we do not include the additional systematic error 
added in (|5.13|) . the standard BBN prediction is found to be more than 2a away from the 
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center value. If we take this discrepancy seriously, we need some explanation which may 
include some effect of particle-physics model beyond the standard model |1H1 EHl HDl HI] • 
Before closing this section, we comment on this issue. 

If the net production of ''Li can be somehow suppressed by the decay of the long-lived 
particle (like the gravitino), the *^Li discrepancy may be solved. In the past, it was discussed 
that the ^Li abundance may be reduced by the photo-dissociation process induced by the 
radiative decay of the long-lived particle |1H] • The scenario with a long-lived particle which 
decays only radiatively is, however, severely constrained by the ^He constraint; in such a 
scenario, photo-dissociation of background *^He is also induced which overproduces ^He. 
Thus, such a scenario does not work once the constraint on the ^He abundance is taken 
into account [301 HI] • 

In order to solve the discrepancy, recently it was pointed out that the suppression of the 
^Li abundance may be possible with hadronically decaying long-lived particles. In |T3] , it 
was discussed that, when the lifetime of the long-lived particle is ~ 10^ sec, abundance of 
'^Li can become consistent with the observational constraint (with no additional systematic 
error) without conflicting other constraints. The reduction of ^Li is mainly due to the 
dissociation of ^Be (which decays into ''Li) by slow neutrons produced in the hadronic 
shower. (To be more exact, such slow neutrons are supplied by the destruction of ^He, the 
inter-conversion from protons, and so on.) 

To study this issue, we have looked for the parameter region where the ^Li abundance 
becomes consistent with the observational constraint 1^5. 12^ ."^^^ Although dissociation of 
^Be by slow neutrons, which is the most important process, is taken into account in our 
numerical calculation, we could not include one of the non-thermal production process 
of ^Li: N + ttBG N + a + tt's, followed by a -|- obg ~^ ''Li + ■ ■ ■ . This is because 
the experimental data for the energy distribution of the final-state a is not available for 
the first process. Thus, the ^Li abundance should be somehow underestimated. With 
a mild assumption that the kinetic energy of the energetic a produced by the process 
+ q;bg —^N + a + it's is ~ 140 MeV in the center of mass system independently of the 
energy of the beam nucleon jSIj, however, we checked that the resultant ''Li abundance is 
not significantly affected by this non-thermal production process in the parameter region 
in which we will be interested. Thus, we believe that our calculation gives a reasonable 
estimate of the ^Li abundance (even though the lower bound on the ''Li abundance was 
not considered in deriving the upper bound on Tr). 

In Fig. ^2 we show the region where the "^Li becomes consistent with the observational 
constraint (j5.12j) . In this calculation, we have used the mSUGRA parameters for the Case 
1 to determine the MSSM parameters although the result is insensitive to the choice of 
the mSUGRA parameters. As one can see, when m3/2 ~ a few TeV and Tr ~ 10^"'' GeV 
(F3/2 ~ 10"^^ - 10-15 ), the ^Li abundance becomes consistent with Eq. ()5.12j) without 

#^^Note in this case that the of logiQ[{mu/n}i)] is calculated by Eq. (not Eq. (|6.2|) '). and 95 

% C.L. corresponds to = 3.84. In addition, because we fix the observational value of nej^i/mj^i, and 
we adopt the observational value of mu/nii in H5.12(l . the observational constraint on nei^^/nu is also 
modified to be {neu/nu)"^^ < {l-WtlH) x 10"" {2a). 
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Figure 11: Parameter region which predicts the ''Li abundance consistent with ()5.12|) : in 
the shaded region the ^Li abundance becomes consistent with ()5.12|1 . We have used the 
mSUGRA parameters for the Case 1. 



conflicting the observational constraints for other hght elements. We have checked that 
this region is consistent with the parameter region suggested in |2] • We have also checked 
that we can flnd a parameter region which predicts ''Li abundance consistent with ()5.11j) . 

7 Conclusions and Discussion 

In this paper, we have studied the effects of unstable gravitino on the BBN in detail. In 
particular, compared to the previous works, we have performed the precise calculations 
of the decay rate and the branching ratios of the gravitino. For this purpose, we have 
flrst flxed the masses and the mixing parameters of the MSSM particles, then calculated 
the decay rates for all the relevant two and three body decay processes of the gravitino. 
Then, we calculate the spectrum of the hadrons (in particular, p, n, and vr^). With the 
hadron spectrum as well as the visible energy emitted from the decay of the gravitino, 
we calculate the light element abundances as functions of the gravitino mass and the 
reheating temperature. By comparing the results of the theoretical calculation with the 
observational constraints, we derived the upper bound on the reheating temperature after 
the inflation. 
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Although we have considered several difference mass spectrum of the MSSM particles, 
the resultant constraints on the reheating temperature behave qualitatively the same. The 
detailed bound is, however, sensitive to the mass spectrum of the superparticles and the 
upper bounds on Tr for several cases are summarized in Table El When the gravitino mass 
is a few TeV, in particular, the hadro-dissociation processes provide significant constraints. 
Of course, in some case, production of the hadrons are suppressed; in particular, when 
the gravitino mass is close to the LSP mass, the only possible two body decay process 
is ip^ — i> 7Xi- In this case, hadrons are produced by the three body decay processes 
ipfi — > QQXiy which is suppressed compared to the two body decay process. Consequently, 
constraints become less stringent. 

If the gravitino is the LSP, the gravitino becomes stable and the cosmological con- 
straints change drastically |21ISI- Detailed study of the case of the gravitino LSP will be 
given elsewhere 
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A Vertex Factors for the Gravitino Decay 



In this appendix, we present the vertex factors for the decay processes of the gravitino. 
For the two-body decay processes with a gauge boson in the final state, the decay rate can 
be calculated with Eq. ((221) with Eq. (jSHI) or (jSSI). The vertex factors &p , Cjf^ c'^\ 

(JT\ 

and Cj^ depends on the mixing parameters. 
For tp^ — > 7x?, the vertex factors are given by 
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(G) 



C 



(G)* 
R 



{92[U*M + gi[U\U 



(A.l) 



where g2 and gi are gauge coupling constants for the SU{2)l and U{1)y gauge group, 
respectively, while for the gluon-gluino final state. 
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and for — > W'^xfy 
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For other processes, the decay rates can be calculated with Eq. ()3.9|) . The vertex factor 
for the neutral Higgs emission processes {ip^ hXi, — > and ip^ — > Ax^) are given 
by 
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respectively. For the rest of the processes (with a quark or a lepton in the final state). 
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B Approximated Formula for Three-Body Process 

Even though there are several diagrams contributing to the three-body decay process of 
the gravitino — > qqxi^ photon-mediated diagram. Fig. [U^a), plays the most important 
role when m^i2 — m^o < mz] indeed, in this case, the decay rate r('?/'^ — > qqxi) is well 
approximated by the results only with Fig. Ql^a). Thus, although we have calculated all 
the relevant diagrams for the three-body processes in our numerical study, we present the 
approximated formula for the differential decay rate for the three-body process, which is 
given by 

dT{^^ ^ qqX^) dVitP^ ^ YXi ^ QQX^) Nc 



qq 9Xi 



2567r3m3/2M2 



X \M\' 



(B.l) 
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where Nc = 3. In addition, m^^ and o are the invariant masses of the qq and gxi 



11 qx'i 

systems, respectively, and are in the following range 



{2mgY < m^- < (m3/2 - m^o)^ 
[rrig + m^oY < mj^o < (^3/2 - mg)^ 



"^3/2 + + 2mJ - (m3/2 - mg) < mj^ + mj^o < ^3/2 + mjo + 2mJ - (m, + m^o) 



(B.2) 

(B.3) 

2 

? 

(B.4) 



with nig being the mass of the final state quark. The photon-mediated three-body decay 
amplitude is given by 
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where p, q', k, k', q are the momenta of ipfj,, Xi, q, q, and intermediate photon 7*, respec- 
tively. Furthermore, C^/^\C^^ are defined in Appendix |Xl and eQg is the electric charge 
of q. 

As discussed in Sectional when 7773/2 — 777^0 > rriz, the process QQXi is mostly 

mediated by the process with on-shell Z-boson, and hence r('?/'^ —>■ qqxi) is well approxi- 
mated by T^ip^ Zxi) X Br(Z qq). 



References 

[1] S. Weinberg, Phys. Rev. Lett. 48, 1303 (1982). 

[2] T. Moroi, H. Murayama and M. Yamaguchi, Phys. Lett. B 303, 289 (1993). 

[3] J. L. Feng, A. Rajaraman and F. Takayama, Phys. Rev. Lett. 91, 011302 (2003); 
Phys. Rev. D 68 (2003) 063504; J. L. Feng, S. f. Su and F. Takayama, P hys. Rev. 
D 70, 063514 (2004); L. Roszkowski and R. Ruiz de Austri, |arXiv:hep-ph/0 408227| 



K. Jedamzik, M. Lemoine and G. Moultaka, arXiv:hep-ph/0506129 



[4] C. L. Bennett et al, Astrophys. J. Suppl. 148, 1 (2003). D. N. Spergel et al, Astro- 
phys. J. Suppl. 148, 175 (2003). 

[5] D. Lindley, Astrophys. J. 294, 1 (1985); M. Y. Khlopov and A. D. Linde, Phys. 
Lett. B 138, 265 (1984); J. R. Ellis, J. E. Kim and D. V. Nanopoulos, Phys. Lett. 



29 



B 145, 181 (1984). R. Juszkiewicz, J. Silk and A. Stebbins, Phys. Lett. B 158, 463 
(1985); J. R. Ellis, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 259, 175 (1985); 
J. Audouze, D. Lindley and J. Silk, Astrophys. J. 293, L53 (1985); D. Lindley, 
Phys. Lett. B 171, 235 (1986); M. Kawasaki and K. Sato, Phys. Lett. B 189, 23 
(1987); R. J. Scherrer and M. S. Turner, Astrophys. J. 331, 19 (1988); J. R. Ellis, 
G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, Nucl. Phys. B 373, 
399 (1992). 

[6] M. Kawasaki and T. Moroi, Prog. Theor. Phys. 93, 879 (1995); Astrophys. J. 452, 
506 (1995). 



[7; 

[8 
[9 
[10 

[11 
[12 

[13 

[14 
[15 
[16 
[17 

[is; 

[19 

[20 
[21 

[22 
[23 



M. Kawasaki and T. Moroi, Phys. Lett. B 346, 27 (1995). 

R. J. Protheroe, T. Stanev and V. S. Berezinsky, Phys. Rev. D 51, 4134 (1995). 
E. Holtmann, M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 60, 023506 (1999). 
K. Jedamzik, Phys. Rev. Lett. 84, 3248 (2000). 

M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 63, 103502 (2001). 
K. Kohri, Phys. Rev. D 64, 043515 (2001). 

R. H. Cyburt, J. R. Elhs, B. D. Fields and K. A. Olive, Phys. Rev. D 67, 103521 
(2003). 

K. Jedamzik, Phys. Rev. D 70, 063524 (2004). 

M. Kawasaki, K. Kohri and T. Moroi, |arXiv:astro-ph/0402490l 

M. Kawasaki, K. Kohri and T. Moroi, Phys. Rev. D 71, 083502 (2005). 

J. R. Ellis, K. A. Olive and E. Vangioni, |arXiv:astro-ph /0503023' 

N. Okada and O. Seto, Phys. Rev. D 71, 023517 (2005). 

K. Kohri, M. Yamaguchi and J. Yokoyama, Phys. Rev. D 70, 043522 (2004); 
arXiv:hep-ph/050221 1, 

M. H. Reno and D. Seckel, Phys. Rev. D 37, 3441 (1988). 

S. Dimopoulos, R. Esmailzadeh, L. J. Hall and G. D. Starkman, Astrophys. J. 330, 
545 (1988); Phys. Rev. Lett. 60, 7 (1988); Nucl. Phys. B 311, 699 (1989). 

H. E. Haber and G. L. Kane, Phys. Rept. 117, 75 (1985). 

P. Gondolo et ai, JCAP 0407, 008 (2004). 



30 



[24] T. Moroi, 'arXiv^hep-ph/9503210l 

[25] The ISAJET homepage: ,http : //mm . phy . bnl . gov/~ isaj et/[ 

[26] PYTHIA webpage: http : / /www . thep . lu . se/~ t orb j orn/Pythia . htiiil[ 

[27] M. Bolz, A. Brandenburg and W. Buchmuller, NucL Phys. B 606, 518 (2001). 

[28] K. Enqvist and M. S. Sloth, Nucl. Phys. B 626, 395 (2002); D. H. Lyth and 
D. Wands, Phys. Lett. B 524, 5 (2002); T. Moroi and T. Takahashi, Phys. Lett. 
B 522, 215 (2001) [Erratum-ibid. B 539, 303 (2002)]; T. Moroi and T. Takahashi, 
Phys. Rev. D 66, 063501 (2002). 

[29] D. Tytler, X. m. Fan and S. Buries, Nature 381, 207 (1996); S. Buries and D. Tytler, 
Astrophys. J. 499, 699 (1998); S. Buries and D. Tytler, Astrophys. J. 507, 732 
(1998); J. M. O'Meara et al., Astrophys. J. 552, 718 (2001); M. Pettini and 
D. V. Bowen, Astrophys. J. 560, 41 (2001); D. Kirkman et al., Astrophys. J. Suppl. 
149, 1 (2003). 

[30] G. Sigl, K. Jedamzik, D. N. Schramm and V. S. Berezinsky, Phys. Rev. D 52, 6682 
(1995). 

[31] J. Geiss, "Origin and Evolution of the Elements," edited by N. Prantzos, E. Vangioni- 
Flam, and M. Casse (Cambridge University Press, 1993) 89. 

[32] B. D. Fields and K. A. Olive, Astrophys. J. 506, 177 (1998). 

[33] Y. I. Izotov and T. X. Thuan, Astrophys. J. 602, 200 (2004). 

[34] K. A. Olive and G. Steigman, Astrophys. J. Suppl. 97, 49 (1995). 

[35] K. A. Olive, E. Skillman and G. Steigman, Astrophys. J. 483, 788 (1997). 

[36] Y. I. Izotov, T. X. Thuan and V. A. Lipovetsky, Astrophys. J. Suppl. 108, 1 (1997). 

[37] K. A. Olive and E. D. Skillman, Astrophys. J. 617, 29 (2004). 

[38] S. G. Ryan, J. Norris, and T. C. Beers, Astrophys. J. 523, 654 (1999); S.G. Ryan et 
al., Astrophys. J. Lett. 530, L57 (2000). 

[39] P. Bonifacio et al., Astron. Astrophys. 390, 91 (2002). 

[40] B. D. Fields, K. Kainulainen, K. A. Olive, and D. Thomas, New Astron. 1, 7 (1996). 

[41] V. V. Smith, D. L. Lambert, and P. E. Nissen, Astrophys. J. 408, 262 (1993); 
L. M. Hobbs and J. A. Thorburn, Astrophys. J. 491, 772 (1997); V. V. Smith, 
D. L. Lambert, and P. E. Nissen, Astrophys. J. 506, 923 (1998); R. Cayrel et al., 
Astron. Astrophys. 343 923 (1999). 



31 



[42] D. K. Duncan, D. L. Lambert, and M. Lemke, Astrophys. J. 401, 584 (1992); 
M. Casse, R. Lehoucq, and E. Vangioni-Flam, Nature 373, 318 (1995); R. Ramaty, 
B. Kozlovsky, and R. E .Lingenfelter, Astrophys. J. Lett. 438, L21 (1995). 

[43] M. Lemoine, D. N. Schramm, J. W. Truran, and C. J. Copi, Astrophys. J. 478, 554 
(1997). 

[44] B. D. Fields and K. A. Ohve, New Astron. 4, 255 (1999). 

[45] J. R. Elhs, G. B. Gelmini, J. L. Lopez, D. V. Nanopoulos and S. Sarkar, in [H]. 

[46] W. Hu and J. Silk, Phys. Rev. Lett. 70, 2661 (1993). 

[47] D. J. Fixsen, et al., Astrophys. J. 473, 576 (1996). 

[48] J. L. Feng, A. Rajaraman and F. Takayama, in |3]. 

[49] K. Ichikawa and M. Kawasaki, Phys. Rev. D 69, 123506 (2004). 

[50] K. Ichikawa, M. Kawasaki and F. Takahashi, Phys. Lett. B 597, 1 (2004). 

[51] C. Y. Wong, "Introduction To High-Energy Heavy Ion Collisions," (World Scientific, 
1994); M. Kaneta [NA44 Collaboration], Prog. Theor. Phys. Suppl. 129, 167 (1997); 
M. Kaneta, Ph. D thesis (Hiroshima University, 1999). 

[52] K. Kohri, T. Moroi and A. Yotsuyanagi, in preparation. 



32 



